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Abstract. The radial distribution of dust has been calcu- 
lated in outflows driven by radiation pressure in which the 
gas drag is taken into account. It is shown that the dust 
distribution in the newly discovered candidate disk around 
BD+31°643 can be explained by this simple model, pro- 
vided that the ambient gas of the surrounding cluster 
IC 348 has a density p ~ 10 4 cm -3 . The dust is prob- 
ably produced within 2300 AU from the central binary 
star. 
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1. Introduction 

Dust around main sequence stars has been discovered 
more than a decade ago by the IRAS satellite (Aumann et 
al. 1984), but little work on the dust dynamics has been 
done. Moreover until recently, /3 Pictoris was the only star 
to capture all the attention, because images showed that 
the dust shell was in fact a disk seen edge-on from the 
Earth (Smith & Terrile 1984). These images have given 
unique information on disk morphology and dust distri- 
bution (Artymowicz et al. 1989, Kalas & Jewitt 1995). 
The dust distribution, the intriguing asymmetries or the 
wrapped geometry have been studied thoroughly (Roques 
et al. 1994, Artymowicz 1995, Lecavelier des Etangs et al. 
1996a, 1996b, Mouillet et al. 1997). But it is still not clear 
if these results apply only to the well-known /3Pictoris disk 
or if they describe a general scheme for which f3 Pictoris 
is an example. 

Recently, the situation changed drastically when Kalas 
& Jewitt (1997) discovered a candidate dust disk around 
BD+31°643, although it is still to be confirmed by ob- 
servational follow-up. As in the (3 Pictoris case, the dust 
lifetime around BD+31°643 is smaller than the age of the 
star. There must therefore be a replenishment mechanism. 
The origin of dust (collision or evaporation of parent bod- 
ies?) must be identified. It is interesting to see if the anal- 
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ysis of the dynamics and the distribution of the dust can 
provide new insight on various aspects of the origin and 
history of dusty disks around main-sequence stars. 

If this dust was in fact coming from the interstellar 
medium and dragged by invisible massive gas, the dynam- 
ics of this disk could be very different from the dynamics 
of the dust around (3 Pictoris. But, if this disk was really 
the second member of the (3 Pictoris family for which we 
have an image, then this would shed more light on the 
origin of the (3 Pictoris-like disks (Lissauer 1997). 

In addition, BD+31°643 is a binary star; this could 
have consequences on possible dust production mecha- 
nisms. At least, if this disk is replenished like the /3Pictoris 
one, this shows that the formation of planctesimals is pos- 
sible around binary stars (Kalas & Jewitt, 1997). The evo- 
lution and stability of planetary systems in this environ- 
ment is starting to be investigated (Brunini 1997, Bcn- 
est 1998). But knowledge on these putative parent bodies 
from the observation of dust have to be improved, and a 
better understanding of the dust dynamics is needed. 

In this paper, after a summary of the characteristics of 
the BD+31°643 candidate disk, where it is recalled that 
the dust is expelled by the radiation pressure (Sect. ||), we 
present theoretical calculation on outflows driven by radi- 
ation pressure in which the gas drag is taken into account 
(Sect. | and |). Application to BD+31°643 (Sect. |) and 
numerical simulations (Sect. |^) show that the observed 
dust distribution (or more precisely the surface brightness 
distribution along the midplane of the disk) can be ex- 
plained by the ambient gas drag of the cluster IC 348, if 
the gas density is between 10 3 and 10 5 cm -3 . 



2. The BD+31°643 disk 

The apparent dust surrounding BD+31°643 presents 
strong similarities with the /3Pictoris circumstellar matter, 
namely a disk-like structure seen edge-on with a scattered 
light distribution close to a power law. The similarities are 
clearly visible through a comparison of the morphology of 
both disks seen on the corresponding images (see Fig. 1 
in Kalas & Jewitt, 1997). 
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Fig. 1. The ratio of the radiation force to the gravita- 
tional force as a function of particle size for "Draine & 
Lee" silicates (dashed line), graphite (solid line) and the 
geometrical optics (dotted line) 

But important differences arc also found. The observed 
dust is at significantly larger distances in the case of 
BD+31°643, simply because the central star is much far- 
ther (~ 330 pc). In the surrounding of BD+31°643, the 
disk is seen between 1300 AU and 6000 AU with a flat- 
ter brightness distribution: F(r) oc r~ 3 between 4 and 
7 arcsec (1300-2300 AU), and F(r) oc r" 19 between 7 
and 20 arcsec (2300-6000 AU). 

Most importantly, the large albedo of the dust in the 
blue shows that this particular dust is smaller than the 
wavelength with a size ~ 0.1/im. As mentioned by Kalas 
& Jewitt, these small grains must be more sensitive to the 
radiation pressure than to the star gravity and must then 
be flowing out from the system. Indeed, the central star 
is a binary system more massive and more luminous than 
/3Pictoris (2x5M Q , 2x83OL ); the ratio of the stellar lu- 
minosity to the stellar mass is then 35 times larger for 
BD+31°643 than for /3Pictoris. Particles of size ~ 0.1/im 
have a (3 ratio of radiation force to gravitational force typi- 
cally 3 to 10 times larger than particles of 1/im size (Burns 
et al., 1979; Artymowicz 1988). If we simply scale the re- 
sults of Burns et al. (1979) for the solar system and of 
Artymowicz (1988) for f3 Pic to the mass-luminosity ratio 
of BD+31°643, we find that particles around BD+31°643 
have (3 ratios between 100 and 350. 

In fact, the (3 ratio could even reach ~ 1000 due to the 
large flux in the UV from a B5V star. With the optical 
properties given by Draine & Lee (1984) for silicates ("DL 
silicates") and for graphite, we carried out the (3 ratio as 



a function of the size of particles (Fig [j]). We have taken 
the UV spectrum measured by IUE on two B5V stars 
used as references (HD 34739 and HD 199081). The radia- 
tion pressure coefficient averaged over the stellar spectrum 
< Qpr > has a major contribution from the flux below 
1500A where the absorption efficiency barely depends on 
the particle size for particles larger than 0.01/im. Conse- 
quently, we find that, in contrast to the solar system case, 
the (3 ratio does not reach its maximum at about 0.1/im 
but around 0.01/im. This shows that particles seen around 
BD+31°643 must have (3 ratio larger than ~100 and must 
be ejected on hyperbolic orbits by radiation pressure. 

Gas is also present within the young cluster IC 348 
in which BD+31°643 is embedded. Its presence is de- 
duced from spectroscopic observations of absorption lines 
and from radio maps (Snow et al. 1994). Snow et al. 
showed that the column density of gas must be Nhi ~ 
3 • 10 21 cm~ 2 , which is also consistent with the extinc- 
tion of the star (£ B -y=0.84, A v = 2.8). This allow to 
obtain observational limits on the gas density in the cir- 
cumstellar environment. With an observed maximal size 
of ~ 1 pc for the IC 348 cluster, the gas density has a 
lower limit pi ~ 10 3 cm~ 3 . The spectroscopic observations 
of Snow et al. show that BD+31°643 is inside the cloud. 
It is also unlikely that the gas density in the star envi- 
ronment is smaller than the mean value in the cluster, 
because BD+31°643 is at the center of the cluster and co- 
incides with the peak of dust extension and gaseous radio 
emission map. On the other hand, if the whole column 
density of the gas was concentrated into the observed disk 
of size > 2000 AU, we would obtain an upper limit of the 
gas density p u ~ 10 5 cm~ 3 . 

In short, the observations show that the disk of 
BD+31°643 has a particular environment with a gaseous 
medium which density is between 10 3 cm -3 and 10 5 cm -3 . 

3. Outflow 

When a dust particle with (3 ^> 1 is produced, it is ejected 
from the system by radiation pressure. It reaches rapidly 
a constant asymptotic velocity Voo ~ vq^/2/3 — 1, where 
vq is the initial velocity of the parent body on its keple- 
rian orbit. Consequently, if some dust is produced close 
to the star, it is spread in a disk-like structure with a ra- 
dial density n oc r~ 2 . Nakano (1990) has given the first 
order relation between the radial distribution of the dust 
density n(r) and F the apparent surface brightness along 
the midplane of the disk if seen edge-on from the Earth: 
F oc n(r)r -1 . As a result, the disk produced by outflowing 
dust has a radial brightness F oc r~ 3 . 

But, as already noted in Lecavelier des Etangs et al. 
(1996a), the gas drag in an edge-on disk decreases the 
gradient of the observed brightness. The next section is 
dedicated to the evaluation of the dust distribution taking 
into account the gas drag. 
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Fig. 2. Plot of the dust radial velocity (dr/dt) as a function of the distance to the star (r) for different gas density 
(p). The dust is supposed to be produced at a distance ro = 100 AU from the star, the ratio of the radiation force 
to the gravitational force is (3 = 200. The full Eq. [I] has been numerically solved (solid line). The fit to the different 
phases: acceleration (dot-dashed lines), deceleration (long-dashed line), stationary outflow (short-dashed and dotted 
lines) are calculated with the approximation given by the equations of Sect ^. 

For p — 30cm~ 3 , the particle is accelerated between 100 and 1000 AU; it is then decelerated between 1000 and 10 5 AU 
and reaches the stationary outflow beyond that point. 

For p — 300cm -3 , the acceleration stops at 300 AU, and the deceleration takes places between 300 and 5000 AU; the 
stationary outflow starts at r-i — 5000 AU. The distribution changes again at r% = 10 5 AU when the thermal drag 
becomes dominant. 

For p = 3000cm -3 , the deceleration phase is very short (between n = 150 AU and r2 = 400 AU). The stationary 
outflow gives a radial velocity f oc r -1 (between r2 = 400 AU and r$ =30 000 AU). The velocity regime changes at 
30 000 AU beyond which we have f oc r -2 . 

For p — 300 000cm -3 , there is no deceleration phase at all. The thermal drag is dominant when r > r% = 2800 AU. 
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4. Outflow with gas drag 

The motion of one particle driven by the radiation pres- 
sure and dragged by the ambient gas is 

dv ln , GM ^ 

^ = (P-V—r + F PR + F g (1) 

where r and v are the position and the velocity of the 
particle, M is the mass of the central star, Fpr and F g are 
the Poynting-Robertson and the gas drags. The Poynting- 
Robertson effect is negligible compared to the gas drag as 
soon as the gas density is larger than ~ 1cm -3 . The gas 
drag is given by (Kwok 1975): 



(a /m)mupp\Jv T + (v - v g ) 2 (v - v g ) 



(2) 



where a and m are the dust cross section and mass, p is the 
mean molecular weight of the gas and ran is the hydrogen 
mass, vt is the sound velocity. v g is the gas velocity and 
is very small, we have taken v g = 0. 

For dust grains, we have a/m = 3/(4sp^) where s 
and p& are the dust size and density. We have shown 
that s oc /3 _1 , where sf3 ~ 20/xm (Sect ||). Then, with 
Pd ~ 1 g cm" 3 , we obtain a/m ~ 375/3 cm 2 g" 1 . 

From the IRAS 60/im emission we know that the 
star environment must be at a very low temperature: 
we can assume T ~ 30 K (Snow et la. 1994). Then, 
as soon as the particle velocity is larger than vt = 
0.14V(T/30K)/^ AU year , the first term under the 
square root of Eq. || is negligible, and assuming p, = 1, 
we get 

F g » -375(3m mP \(v - v g )\(v - v g ) » -C/3p\v\v (3) 
where C = 9.39 • 10" 9 cm 3 AU" 1 . 

In order to evaluate the radial distribution of dust, we 
carry out the different phase of the motion by approximate 
solutions of Eq. |l|. We consider only the radial component 
in the polar coordinates (r, 6) of the vectorial Eq. pj: 



<? r , de ,2 ro ^ GM 



(4) 



Since we know that n(r) ~ f _1 r -2 , we simply need to 
determine the radial velocity r as a function of the distance 
to the star, r. Then, the apparent surface brightness along 
the midplane of the disk, F oc n(r)r -1 , can be compared 
to the observations. If (3 3> 1, the Eq. || becomes 

d 2 r 



dt 2 



(3 



GM 



(5) 



4-1- Acceleration phase 



When a particle is produced, its radial distance and ve- 
locity are small. Consequently, the last term of Eq. ^| is 
negligible and we obtain 



d 2 ^ 
dt 2 



GM 



This gives f — y / 2f3GM y r o r 1 , where ro is the ini- 
tial radial distance at which the dust was produced. This 
approximation is plotted on Fig. ^ with dot-dashed lines. 
This is valid for small velocity and small distances, as long 
as d 2 r/dt 2 > 0, that is to say for r < n where 



n = (Vo + ^r 2 + 2r /Cf3p^j /2. 



(7) 



.2. Deceleration phase 



After the particle has gone farther from the star, where the 
radiation pressure becomes less important, and has gained 
a large radial velocity, the gas drag becomes dominant and 
we have 



dt 2 



(8) 



During this deceleration phase, the first term within 
parenthesis is largely dominant, and the second term can 
be considered to be a constant where r 2 f 2 is approx- 
imated by its mean value < r 2 r 2 >. Hence, we have 



-ix{r—r\) 



This approximation is plotted on Fig. ^ with long- 
dashed lines. This is valid as long as the velocity is large 
enough for the drag to dominate the radiation pressure. 
Subsequently, the radial velocity significantly decreases 
and the drag decreases to the radiation pressure's level. 
Thus the deceleration takes place when r\ < r < r2 , where 
r2 is defined by 



-a(r 2 -r 1 ) 



I GM 

-op- 1 



(9) 



4-3. Stationary outflow 

When r > r2, then we simply have d 2 r/dt 2 ~ 0: this is 
the stationary outflow where the gas drag exactly com- 
pensates the radiation pressure. 

From Eq. ||, we get r = y^jjf f 1 - This gives 

n(r) oc and F(r) oc r~ 2 , exactly as observed around 
BD+31°643. 

This is represented on Fig. ||by the short-dashed lines. 



When the particle velocity becomes very small, the 



thermal term is dominant in the gas drag (Eq. 
situation occurs when r > r 3 where r 3 = 
is then 




GM 



This 
Eq.| 

(10) 



where C' = 375to h AWt ~ 1-3 • 10- 9 v /(T/30K)cm 3 year- 1 . 
During this phase, we still have d 2 r/dt 2 ~ 0, which gives 
2 , and F{r) oc r _1 . This approximation is plot- 



- GM , 
C'p 



(6) ted on Fig. || with dotted lines. 



A. Lecavelier des Etangs et al.: Dust distribution in the BD+31°643 disk 



5 



: o 4 



1000 



100 




10 



r. = 10 jS=100 



-10 /S-200 
: 1 /? = 200 
1000 ^-200 



r =1000 /3=400 



100 



1000 



10 4 10 5 
p(crrr 3 ) 



10 6 



Fig. 3. Plot of the radial distance r 2 , as denned in the text. 
The hashed zone represents the region where r > r 3 , for 
a temperature T = 30 K. Above the lines representing r 2 



and below the hashed zone, we have t 



hence F(f) oc 



r as observed around BD+31°643 between 2300 and 
6000 AU. We see that r 2 < 2300 AU gives p > 10 3 cm-\ 
and r 3 > 6000 AU gives p < 10 5 cm _1 , 



5. Application 

5.1. The radial velocity versus gas density 

We applied the previous calculations and approximations 
to the case where a grain is produced at a distance 
r = 100 AU from the star and with j3 = 200. We have 
plotted the dust radial velocity (dr/dt) as a function of 
the distance to the star (r) for different gas density (p) 
(Fig. ||). The full Eq. [I] has been numerically solved, and 
the result is given by the solid line in Fig. |^. 

The fit to the different phases: acceleration (dot- 
dashed lines), deceleration (long-dashed line), stationary 
outflow (short-dashed and dotted lines) are calculated 
with the approximation given by the equations of Sect ^. 
The superposition of the exact solution of Eq. [l] and the 
different fits validates these approximations and the cor- 
responding domain of application given in the previous 
section. This also validates the different positions where 
the changes of regime occur (r%, r 2 , 7-3). 

For p = 30cm~ 3 , we see that the particle is accelerated 
between 100 and 1000 AU; it is then decelerated between 
1000 and 10 5 AU and reaches the stationary outflow be- 
yond that point. 

For p = 300cm -3 , the acceleration stops at 300 AU, 
and the deceleration takes places between 300 and 
5000 AU; the stationary outflow starts at r 2 — 5000 AU. 



The distribution changes again at r 3 = 10 5 AU when the 
thermal drag becomes dominant. This gas density is too 
low to explain the observed distribution of dust around 
BD+31°643. 

For p = 3000cm -3 , the deceleration phase is very short 
(between r x = 150 AU and r 2 = 400 AU). The sta- 
tionary outflow between 400 AU and (r 3 =)30 000 AU 
gives a radial velocity f oc r . Then the disk surface 



brightness distribution is F(r) 



, exactly as ob- 



served around BD+31°643. The distribution changes at 
30 000 AU beyond which we have f oc r~ 2 , and conse- 
quently F(r) oc r _1 . 

For p = 300 000cm -3 , there is no deceleration 
phase at all. The thermal drag is dominant when 
r > r 3 = 2800 AU. This gas density is too large to explain 
the observed distribution of dust around BD+31°643. 



5.2. The BD+3F643 disk 

In the BD+31°643 candidate disk, the distribution ob- 
served between 2300 and 6000 AU can be explained if 
the outflow is stationary between r 2 and r 3 as described 
in Sect. 4.3. This happens only if r 2 < 2300 AU and 
r 3 > 6000 AU. In Fig. || we have plotted these dis- 
tances for different ratio f3 and initial distance rg. We 
can see that these conditions are satisfied if we have 
10 3 cm -1 < p < 10 5 cm -1 . 

Incidentally, another solution could be found if the 
dust is in the stationary outflow in the thermal drag 
regime (r > r 3 ) and with a radially decreasing gas density 
p oc r _1 . In that case, we have r oc r~ 2 j p oc r" 1 , and 
again F oc r~ 2 . But this would imply that the gas den- 
sity should be larger than 10 6 cm -3 which is incompatible 
with the observed gaseous column density. 



6. Numerical results 

For a mixture of particle size (and f3), a radial distribution 
of the sources and the exact treatment of the equation of 
motion (Eq. |l|), numerical simulations are required. We 
built a code to solve this equation and to calculate the 
distribution of the dust particles in the steady state of 
the equilibrium between the continuous replenishment by 
a distribution of parent bodies and the ejection by the 
radiation pressure. 

A wide range of parameters for the dust size distribu- 
tion, the gaseous density and velocity (v g ), and the initial 
parent bodies positions have been used. The simulations 
show that dust produced within 2300 AU is able to ex- 
plain the surface brightness of the disk, provided that the 
gas density around the star is p ~ 10 4 cm -3 . An exam- 
ple of such a simulation is given in Fig. 0, with a size 
distribution corresponding to (3 S [10, 1000]. Any size dis- 
tribution gives the same results, if particles with /3 10 
are in negligible quantity. 
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Fig. 4. Brightness distribution of a dust disk seen edge-on 
and produced by parent bodies on circular orbits between 
2000 and 2300 AU. The gas drag is taken into account with 
p = 10 4 cm~ 3 . The particles have a (3 ratio of the gravi- 
tational forces to radiation forces between 10 and 1000. 
This distribution is very close to the surface brightness 
distribution observed around BD+31°643 by Kalas & Je- 
witt (1997) (dashed line, slightly arbitrarily shifted to be 
better compared.) 



All these parameters are very consistent with the ob- 
servations of the large albedo of grains in the blue and 
the estimated density of gas at the center of the cluster 
IC 348 (Sect |2|). The only parameter newly constrained 
by these simulations is the position of the parent bodies, 
which must be located between 2000 and 2300 AU from 
the star to explain the observed decrease of slope within 
2300 AU. It must be mentioned here that BD+31°643 is 
200 times brighter than j3 Pictoris, so that the distance of 
evaporation of CO (150 AU around /3Pictoris) corresponds 
to about 2100 AU around BD+31°643. But this change of 
slope can also be linked to the limited resolution of the 
observations and needs confirmation. 

It must also be noted that the velocity of the star rela- 
tive to the ambient gas (i>* ) could modify the appearance 
of the disk. However, as can be inferred from the Fig. pL 
simulations show that asymmetries would be visible only 
for distances r k, 10000 (v*/lkm s _1 ) _1 AU. The obser- 
vations are thus consistent with the present modeling if 
u* Si 2km s — . Considering that the star is at the center of 
the cluster and that only the transverse motion would give 
noticeable asymmetries, this last condition is certainly not 
unlikely. 



Finally, the gas drag is the dominant phenomenon 
which allows to explain the dust radial distribution, but 
this also rules out the possibility to better constrain the 
dust origin. 

7. Discussion 

The dynamics of particles in the BD+31°643 disk and 
the observed distribution of dust can be explained by the 
simple model of dust outflow presented here. 

Unfortunately, the issue of the origin of the very small 
particle size remains unexplained. It is not clear why we 
see so much small particles around BD+31°643, whereas 
the dust around other main sequence stars with known 
infrared excess is usually much larger (e.g., Habing et al. 
1996). 

Also the process which produces these grains is still not 
determined. The present model of the BD+31°643 disk 
does not discriminate between production by collision or 
by evaporation. However, a change of slope and possible 
inner hole are observed at a distance of 2000 AU, corre- 
sponding to the sublimation distance of CO (~ 2100 AU). 
Production of dust by evaporation could start if some bod- 
ies were heated by the central star after the dissipation of 
the dense optically thick disk from which a planetary sys- 
tem formed. Alternatively, parent bodies could evaporate 
if they were scattered planetesimals on eccentric orbits in 
a forming Oort-like cloud or in a Kuiper-like outer belt 
(Scattered Icy Objects? Duncan & Levison 1997). 

Finally, the particles around BD+31°643 can have 
a similar origin and nature as the dust grains around 
(3 Pictoris, but both disks present different characteris- 
tics. The main differences are related to the large dis- 
tance ratio compensated in part by the larger luminosity 
of BD+31°643, and consequently the dynamics of dust 
here expelled by the strong radiation pressure. 
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